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Abstract. We outline a gauge-invariant framework to calculate cosmological perturbations 
in dark energy models consisting of a scalar field interacting with dark matter via energy 
and momentum exchanges. Focusing on three well-known models of quintessence and three 
common types of dark sector interactions, we calculate the matter and dark energy power 
spectra as well as the Integrated Sachs- Wolfe (ISW) effect in these models. We show how 
the presence of dark sector interactions can produce a large-scale enhancement in the matter 
y—i power spectrum and a boost in the low multipoles of the cosmic microwave background 

anisotropies. Nevertheless, we find these enhancements to be much more subtle than those 
found by previous authors who model dark energy using simple ansatz for the equation of 
state. We also address issues of instabilities and emphasise the importance of momentum 
exchanges in the dark sector. 
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1 Introduction 

Observation of type la supernovae at high redshift, large scale structure surveys and ob- 
servations of the cosmic microwave background (CMB) all indicate that the Universe is 
currently undergoing accelerated expansion [1-4]. The most convincing explanation for this 
phenomenon is that the Universe is filled in great abundance with a form of energy with neg- 
ative pressure. This so-called 'dark energy' (DE) can lead to an accelerated cosmic expansion 
if its equation of state w satisfies w < —1/3, although observations currently constrain w 
to be close to —1. At present, the cosmological constant model in which w is exactly — 1 
appears to be the most robust model for dark energy [2, 5]. 

However, there remain unsolved problems with the physical origin and the exact value 
of the cosmological constant. It is well known that the value of the cosmological constant, 
interpreted as vacuum energy arising from particle-antiparticle annihilation, and that con- 
sistent with observation can differ by as much as 120 orders of magnitude [6]. Moreover, the 
value of the cosmological constant, however miniscule, must be set at its precise value at the 
initial inflationary epoch, with the slightest deviation from this value leading to disastrous 
consequences [7]. 
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Over a decade since the discovery of cosmic acceleration, there are now an overwhelming 
number of dark energy alternatives to the cosmological constant (see e.g. [8] for a review). 
One of the earliest and the simplest of such models is that of a single light scalar-field model 
of dynamical dark energy dubbed 'quintessence' [9—11]. Unlike the cosmological constant, 
many quintessence models exhibit an attractor behaviour and thus do not require initial 
conditions that are extremely fine-tuned. The late-time acceleration can also be dynamically 
'switched on' without requiring a cosmic coincidence [12]. 

Since quintessence can be interpreted as a field with non-zero mass, it is possible that 
dark energy could cluster like dark matter, or even interact with dark matter. The possibility 
that there may be interactions in the dark sector has previously been explored by several au- 
thors [13-28]. It is common in the previous works to either neglect dark energy perturbations 
or model dark energy using some phenomenological form of uu (e.g. w = constant). These 
simplistic approaches may not accurately reproduce the dynamics of quintessence models as 
suggested by [29, 30]. 

The primary aim of this paper is to investigate the prospects of detecting the observa- 
tional signatures of dark sector interactions. We shall consider three well-known models of 
quintessence and study the effects of dark sector interactions in both the background and 
in the linear perturbations. In particular, we shall examine the effects of non-zero interac- 
tions on the linear matter power spectrum and the integrated Sachs- Wolfe (ISW) using the 
gauge-invariant approaches of Kodama and Sasaki [31]. Unlike previous works on this sub- 
ject, our approach takes into account dark energy perturbations, including the perturbative 
effects of energy and momentum transfer in the dark sector. Issues of instabilities previously 
discovered in [14, 15, 27] will also be addressed. 



2 Quintessence models 

Quintessence refers to a scalar field, <j), evolving along a potential, V(cp), which becomes 
sufficiently flat at late times, leading to cosmic acceleration (see e.g. [8]). At early times, the 
dark energy contribution to the overall energy density is small [32]. Some classes of potentials 
exhibit tracking behaviour whereby the dark energy density tracks below that of matter and 
radiation and only comes to dominate the universe at late times [33]. Another interesting 
feature of quintessence is the existence of attractor solutions [32] which greatly reduce the 
need for the initial values of the dark energy parameters to be finely tuned. 

We will focus our analysis on three models of quintessence, namely, the Ratra-Peebles, 
SUGRA and double exponential potentials summarised below. From this point on, all equa- 
tions will be in natural units with c = H = 1. We also write R 2 = 8irG. 



2.1 Ratra-Peebles 

Ratra and Peebles (RP) [9] showed that the class of inverse power law potentials of the form 

= ~^T' a > °i 

tracks the equation of state during matter and radiation-dominated eras, and dark energy 
only becomes dominant at late times. This behaviour holds for a wide range of initial condi- 
tions. We will use M = 4.1 x 10~ 28 mpi and a = 0.5 as indicated by [34] who obtained these 
values from a likelihood analysis of the WMAP 5-year data [2] and the 'union' supernovae 
data [3]. 
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2.2 SUGRA 

The SUGRA model with potential 

, s M 4+a exp(^ 2 2 ) , 
V{4>) = p 2 V \ (2-2) 

is motivated by supersymmetry and differs from the Ratra-Peebles (RP) potential by an 
exponential 'supergravity' correction factor [35]. This factor can be shown to push the 
equation of state closer to —1 at late times. We will use the values M = 1.7 x lCP 25 mpi and 
a = 1 as indicated by [34] . Again, the SUGRA potential has attractor solutions and is viable 
for a wide range of initial conditions. 

2.3 Double exponential 

The double exponential (DExp) potential 

V{(j)) = M ie - Xl<t> + M 2 e- X2 ^, (2.3) 

allows dark energy to track the equation of state of radiation and matter and leads to acceler- 
ation at late times [36]. We shall use M x = 10" 14 m P i, M 2 = 10~ 13 m P i, Ai = 9.43 and A 2 = 1 
which, as shown in [37], should give rise to observables that are significantly different from 
the cosmological constant but are still broadly consistent with observational constraints. 

3 Evolution of background energy densities 

We assume an isotropic, homogeneous and spatially flat background as described by the 
Friedmann-Robertson- Walker (FRW) metric, 

ds 2 = -dt 2 + a{f) 2 dx i dx i , (3.1) 

where a(t) is the scale factor as a function of cosmic time, t, and summation is implied 
over i = 1, 2, 3. We assume that the background energy density comprises matter, radiation 
and dark energy (denoted by subcripts m, r and (j) respectively). The total energy density, 
P = Pm + Pr + P4>-> satisfies the Friedmann and acceleration equations 

/ • \ 2 -2 

H2 ={1) =^{prn + Pr + Pc l >), (3.2) 
R 2 ( 4 \ 

H = [p m + ^Pr + P^> + Ptf, j , (3.3) 

as well as the energy conservation equation 

p + 3H(p + p) = 0. (3.4) 

Here overdots are derivatives with respect to t and p x is the pressure of component x related 
to its energy density via the equation of state w x = p x /p x - In the above, we have used 
w m = 0, w r = 1/3 and defined below. 

The quintessence field evolves via the Klein-Gordan equation 

4> + 3Hj> + V'{4>) = 0. (3.5) 
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Its energy density, pressure and equation of state are given by 

= = ? ' 2 - V M. (3.7) 

The case in which = reduces to the cosmological constant with w\ = — 1. Note that at 
present uu is observationally constrained at low redshift to a value close to = — 1 with 
~ 10% accuracy [38, 39]. 

We model interactions between dark matter and dark energy by including a source 
term, Q x , in the evolution equations for the dark energy and dark matter background energy 
densities. 

p x = -3Hp x (l + w x ) + Q x . (3.8) 

Non-zero Q thus represents the energy flow from one component to another. To conserve the 
total energy, we require 

Y^Qx = 0. (3.9) 

X 

If interactions are purely between the dark matter and dark energy (i.e. interactions only in 
the dark sector) then 

Qm = -Q<t>- (3-10) 

3.1 Evolution equations 

To evolve the background energy densities with time, we use the set of dimensionless 'energy' 



variables defined in [33] , 



R V 



vW y #V3' 

r = H\J ~3' m= i?V^' A = "^#- (3 - 11} 

Intuitively, x, y, r and m can be thought of as (the square-root of) the ratio of energy density 

of each component to the total energy density (~ 3H 2 ). Converting the time variable in 
Equations 3.3, 3.4 and 3.5 to the 'e-fold' number, N = In a, we obtain 

dhxH 3 ( , 9 r 2 \ , 

'l + x 2 -y 2 + - , (3.12) 



dN 2 V 3 
dx /3 , 2 ( - d\nH\ k 2 



dN \2 Xy - X { 3 + ^N-) + 6lp-x Q ^ (3 ' 13) 

d ^ = -\f*\xy-y d]n ^, (3.14) 

dN V 2 y y dN ' 1 ; 

dr ( dkiH\ 

2 + -7^ . (3-15) 



dm / 3 dlnH\ K 2 
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Figure 1. Left: The equation of state w ( j,(z) of dark energy for the cosmological constant (CC), Ratra- 
Peebles (RP), SUGRA and double exponential (DExp) models (see section 2). Right: The evolution 
of the density parameter £l x = p x / p for the CC and DExp models. Apparently large deviations from 
w<p = — 1 lead to subtle differences in fi at late times. 

These variables obey the Friedmann constraint derived from dividing Equation 3.2 by H. 

x 2 + y 2 + r 2 + m 2 = 1, (3.17) 

We solved these coupled differential equations by integrating them numerically from 
redshift z = 10 8 until z = (today). Initial values were chosen by requiring that the final 
(z = 0) values of the component energy densities agreed broadly with those measured today: 
n m = 0.27, n r = 8.6 x 1CT 5 , n<j, = 1 - n m - U r and H = 70.8 kms _1 Mpc _1 , where 
&x = Px/p- 

The evolution of for the four models and the background energy densities for the 
CC and DExp models are shown in Figure 1, in which the eras of radiation, matter and 
dark energy domination are clearly seen in that order. The quintessence and CC models all 
satisfy ~ — 1 until z ~ 1. Despite an apparently large deviation in from —1 at very 
low redshift, the differences in the background density between the DExp and CC models are 
surprisingly small. Thus, it may be more fruitful to try to distinguish between quintessence 
and the cosmological constant via perturbations. 

4 Perturbations 

We consider only scalar perturbations to the FRW metric. Working in the Newtonian gauge 
and ignoring anisotropic stresses, the perturbed metric takes the form 

ds 2 = -(1 + 2<5>)dt 2 + a(t) 2 (l - 2<S>)dx i dx i , (4.1) 

where $ is the Newtonian potential (for details see [40]). The components of the energy- 
momentum tensor, are 




T ° = -p, T l ° = --(p + p)v, l , T^=p5/, (4.2) 

where k is the Fourier wavenumber, Vi is the peculiar velocity of the perturbation for species 
i and a comma denotes a partial derivative. Each energy density can be decomposed as 
p = p + 5p where barred quantities are the background values. We then define the variables: 
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d m =^— , O r = — , d(p = <p—(p. 
Pm Pr 

Physically, 5 is a measure of density fluctuations (overdensities and underdensities). To avoid 
unphysical gauge artefacts, we shall consider the gauge- invariant overdensities A x defined as 
[31] ' 



aH 



(1 - q m )v m , (4.3) 



A r = 5 r + A {^Yj (! - Qr)v r , (4.4) 

#^» + (3g0 + y / (0))^-0 2 a> 

A^ = , (4.5) 

P4> 



n a H ^ \ 3,, . / aH\ 2 -r- 
+ 3-(l + ^ +2 (l + .,)(-) A, 



(4.6) 

where t>^ = ^^, (4.7) 
a</> 



3i7p x (l + Wj) ' 



(4.8) 



It will also be necessary to consider the perturbative effects arising from the dark sector 
interactions. We define e to be the perturbation in the energy transfer Q x — > Q x 

Q x = Q x (l + e x ), (4.9) 

and note that e x obeys the conservation equation 

J2Q* e * = °- ( 4 - 10 ) 

X 

Similarly, the gauge-invariant variable for energy transfer, E x , can be constructed as [31] 

E x = e x - °§^-v x . (4.11) 
4.1 Covariant formulation of dark sector interactions 

We now describe how dark sector interactions can be covariantly formulated and incorporated 
into the gauge-invariant system (see [14, 15, 23] for previous attempts in this direction for 
simple dark energy models). 

Dark sector interactions can be covariantly described by the conservation equation 

v„rr = Q%- (4.12) 

The 4- vector Q x can be constructed using the function Q x by 

Q(x)fi = Q(x) u (a)fi, (4-13) 

where u^ a ^ is the 4- velocity of component a (which may be different from x). By perturbing 
(4.13), we find 

Q{x)n = Q(x)Ufi + /(x) M > ( 4 ' 14 ) 
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where is the average velocity. The perturbed average 4-velocity — > and perturbed 
4-velocity of component a (uu\^ — > u^^) can be expressed as 



= u {a)0 = -a{\ + $), 



Uj = avYj, 



av a Yj. 



(4.15) 



where Yj is a basis vector constructed from harmonic functions 1 . In 4.14, the vector f( x )iJt 
can be interpreted as the momentum exchange with components given by 



f(x)o = 0, f( x )j = aHp x (l + w x )f x Y 3 



31 



(4.16) 



where f x is an ordinary scalar representing the amplitude of the momentum exchange. In 
addition, f x satisfies the conservation of momentum 



(4.17) 



As with the energy exchange, we can similarly recast f x in a gauge-invariant form F x 

F, = U- £ xi H~ V) y (4-18) 
Hp x (l + w x ) 

All the ingredients introduced so far are now sufficient to allow us to study the evolution 
of the perturbation variables {A m , A r , As, v m , v r , v^}. In the regime where perturbations are 
small, we find the following set of coupled differential equations [31] 

dA m 



where 



dN 

dA r 

~d~N 

dA 6 
dN 



dv r , 



dN 
dv r 

dN 



dN 



9 ClH n QmA m 



k Qrn,H m 

aH p m H 



all , . . . 
6—(l + w)(v - v r ) + A r 



4k 

3aH^ 



3w 4 

2x*_ 

P<t> 



A A + 



Hp4,J 

9 aH .._ . k 

- T (l + ^--^)--^ 



P<t>H 



+ (1 + w&)F& 



3aH - „ 
-v m - -r^A + F, 



-V r + 



-V,/, + 



2k 

k 

\~a~H L 



(1 + ws)aH 



m j 
3aH 
2k 

A^ 



3aH 
~2~k 



A + Fs 



1 



1 + w 



2x 2 v ( f > + ^r 2 v r + m 2 v m , . 



A = (m 2 A m + r 2 A r + (x 2 + y 2 )A 6 ) + l^Q. 



W = -{p r W r + PAWS) 
P 



+ (x 2 + y 2 )ws- 



(4.19) 



(4.20) 



(4.21) 
(4.22) 

(4.23) 

(4.24) 

(4.25) 
(4.26) 

(4.27) 



see Appendix C of [31] for detail. We will not require the explicit form of Yj in the evolution equations. 
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Note that by setting Q X ,E X and F x to 0, we recover the non-interacting case. The 
initial values of the perturbations are fixed using adiabatic initial conditions [41, 42] which 
we calculate to be 



A, 



1 - q m 



3 A 

4 A r , 



A,. 



A, q x 



Q, 



SaHp x 



k A m 
3aH 1 — q„ 



kA r 



V<f, 



0. 



(4.28) 



The only free variable is the initial matter density perturbation which is determined by 
normalising the matter power spectrum at a pivot scale which we take to be k = O.OOlMpc -1 . 
The power spectrum of energy component x is calculated by comparing the amplitudes of 
the growing mode at early and late times 



Px(k) = k r ' 



A x (z = 0) 



1 2 



where n s is the primordial scalar spectral index taken to be 0.96. 
4.2 Power spectra without interactions 



(4.29) 




Figure 2. Left: If there is no interaction in the dark sector, the linear matter power spectra for all 
four 'best-fit' dark energy models are essentially indistinguishable. Right: the fractional percentage 
difference between the matter power spectra of the quintessence models and the cosmological constant. 
The best-fit models all exhibit sub-percent differences from the cosmological constant. 

Figure 2 shows the matter power spectrum for all four dark energy models as well as 
the fractional difference in the matter power spectra for the quintessence models relative to 
the cosmological constant expressed as a percentage. 

p /L\modcl p (h.\CC 



P m (k) 



cc 



Note that our calculations are not valid for scales k > 0.1 where we expect a significant 
enhancement in the power spectra due to non-linearities. The equations governing the per- 
turbations in the non- linear regime are beyond the scope of our work. To study the perturba- 
tions in this regime, one requires an N— body simulation as in [34, 43] or a more sophisticated 
perturbation theory [44]. 
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Figure 3. The linear power spectra for dark energy for the three quintessence models with no 
interactions. Note the 'turnover' on very large scales close to the Hubble radius (~ 2 x 10~ 4 Mpc -1 ) 
indicating the typical size of dark energy perturbations. 



Differences between the models are small with the largest difference ~ 1% observed for 
the DExp potential. The quintessence models have an excess power for k > O.OOlMpc- 1 and 
a power deficit for k < O.OOlMpc^ 1 . This correlates well with the 'turnover' to the dark 
energy spectra shown in Figure 3. Nevertheless, these small differences of order 1% in matter 
power spectra for the models are essentially unobservable. 

The power spectrum for dark energy perturbations, shown in Figure 3, are roughly four 
orders of magnitude smaller than the matter power spectra. Thus, we do not expect to see 
large differences in the matter power spectra of the quintessence models due to dark energy 
clustering. We see that the large-scale clustering of dark energy enhances the clustering of 
dark matter on large scales. The dark energy power spectra are approximately constant on 
very large scales (k < O.OOlMpc -1 ), and decay exponentially on small scales. In fact, dark 
energy clustering is significant mainly on large scales in all models in which the sound speed 
of dark energy equals the speed of light [45] . The decay of the power spectrum is split into 
two regions with different exponents changing near the turnover scale in the matter power 
spectrum. The differences in the power spectra of the different quintessence models is due 
to the different values of at late times, since A m oc (1 + w) (see 4.19 and 4.27). We will 
shortly investigate the effects of dark sector interactions on these power spectra. 



5 Integrated Sachs-Wolfe effect 

The Integrated Sachs- Wolfe effect (ISW) is an effect observed in the CMB for small multipole 
number t caused by the blueshifting of photons as the gravitational potential wells become 
shallower due to the late-time cosmic acceleration [46]. The gravitational potential $ can be 
decomposed as: 

$ = $ m + $ r + $ . (5.1) 
The potential <I> is related to the perturbation variables by the cosmological Poisson's equation 

3 faH\ 2 



* = "(^rj A - (5 - 2) 
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Figure 4. The percentage differences between Cj for the three quintessence models and the CC 
model. The contributions from quintessence are overwhelmed by cosmic variance, hence the models 
are observationally indistinguishable via the ISW effect. 



Consider the CMB temperature anisotropy AT/T(n) along the line-of-sight unit vector 
n. The CMB anisotropy power spectrum is given by 



/AT AT 
( — (n) — (m) 



nm=[i 



i^+l)Q%), 



(5.3) 



where Ci is amplitude of the CMB angular power spectrum at the £-th multipole and Pe is 
the Legendre polynomial of order I. By decomposing AT/T into the spherical harmonics 
and integrating over all directions we find [47] 



C, 



ISW 



oo r 



k 2 



dec 



a dz 



-Jf k 



1 



H(z')(l 



dz' 



dz dk, 



(5.4) 



where Zdec is the redshift at which photon decoupling occurs (~ 1020) and Je(x) is the 
spherical Bessel function of order I. 

Figure 4 shows the squared percentage differences between Cj for the three quintessence 
models and the cosmological constant model 



(£,quint _ gCC * 

in comparison with the cosmic variance 



x 100%, 



2^ + 1 



(5.5) 



The relative differences in the ISW effect for quintessence models and ACDM are com- 
pletely overwhelmed by cosmic variance. They are also much smaller than indicated in other 
investigations in which the equation of state for dark energy were parameterised by simple 
ansatz e.g. = wq or = wq + w\z (where wq and u>i are constants). This is perhaps 
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not surprising since our 'best-fit' quintessence models can closely replicate ACDM dynamics 
whereas these parametrizations can diverge from ACDM dynamics substantially. 

In summary, we find that the ISW is an ineffective discriminant between the cosmological 
constant and models of quintessence that deviate from w = — 1 at very low redshifts. In 
the next section, we investigate whether this conclusion is changed if there are interactions 
between quintessence and dark matter. 

6 Modelling dark sector interactions 

In this section, we introduce three models of dark sector interactions and explain how they 
can be rewritten in covariant forms. 

6.1 Dark matter decay 

We first consider a dark sector interaction of the form 

Q m = -Ap m , (6.1) 

which represents the decay of dark matter into dark energy with a constant decay rate. The 
decay rate per unit energy density of dark matter is quantified by the positive constant A, 
which we refer to as the interaction strength. This form of interaction has been previously 
investigated in [14]. 

In the covariant formulation, the interaction term is given by Ap m u m , where u m = 
(—a, 0), since this is the rate at which 4- momentum is transferred to dark energy. Hence, 

Q(m)n = -ApmUirn)^- (6-2) 
We perturb the above equation to find 

Q(m)fi = Qm(l + £m)«(m)/i) Q<t> = + ^m) M (m)/i- ( 6 -3) 

Using Equations 4.3, 4.9, 4.11, 4.15 and the above, we find the gauge-invariant energy transfer 
e m = Sm, E m = A m , E^ = A m + 3^ ( 1 - ® m ) (y^ - v m ). (6.4) 



k V 3H P 

Similarly, using Equations 4.14, 4.16, 4.17 and 4.18, we find the momentum transfer 

, _ QmiVm ~ V) p _ p _ Q^Vm-V^)- (R 
Jm - Tj , r m — U, t^- — — -. (b.5j 

6.2 Dark energy decay 

The second model represents the case in which dark energy decays into dark matter at a 
constant decay rate [48]. This may be applicable, for instance, in the early dark energy 
scenario [49]. The interaction term has the form 

Q<j> = -Ap<p, Q m = Ap^, Q(4>)^ = -Ap^u^^. (6.6) 



Similarly, using Equations 4.9-4.18, we find the gauge-invariant energy and momentum trans- 
fers 

+ . (6.7) 



ciH 

E (j> = A -3— (l + «v) 



3^(1 + W<j> ) 2k 
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E m = A (f) + 3 a ^-(l + w 4 



(6.« 



Fm = 9^± <l . (6 .9) 

Note that the momentum transfer F m and Fa, are 'orthogonal' to those in the dark matter 
decay model. 

6.3 Scalar-tensor model 

Finally, let us consider a interaction term which appears in some scalar-tensor theories of 
gravity [50] 

Q{4>)H = ApmV^, Q( m )ti = -ApmV^cj). (6.10) 

Using Equation 4.14 we find 

Qa> = -WAp m {\ + 5 m )d^j), = Ap m 4> (6.11) 

where u^ 1 is the average velocity. To calculate the gauge-invariant interaction, we use the 
following useful expressions for the perturbed scalar field 

& = i-*i + 8l f* = ^, (6.12) 

^ = _N_ + g_ W + VW (613) 

Using Equations 3.9, 4.15 and 4.10, we find the gauge-invariant energy transfer variables 



A<4 . 3aH ( Q m „\ 3 / aH\ 



2 



E * = A " + TT^ + — l^-V-^l-J A ' (6 ' 14) 

T-i a A , 3aH { Q m v m \ 3/ai7\ 2 x 3aH y 2 X 

Similarly, from 4.7, 4.16, 4.18 and 6.10, we find the gauge-invariant momentum transfer 
variables to be 



Ap^^ - v) 

U = ~e- — rr, — v F 't> = °i F m = — f7 • 6.16 

Hp^l + u^) Hp m 

7 Results 

We are now ready to calculate the observables for the different interaction models. In particu- 
lar, we would like to establish if dark sector interactions could give rise to observable imprints 
in the power spectrum or the ISW effect whilst satisfying the constraint that w is very close 
to —1. We use the SUGRA potential for the dark energy decay and scalar-tensor model 
because, at the background level, it behaves most similarly to the cosmological constant (see 
Figures 1 and 2). It is computationally more challenging to implement the SUGRA potential 
with the dark matter decay model due to the non-monotonic nature of the potential and the 
RP potential was used instead (where an analytic solution for <j> can be calculated given the 
value of V((p)). In all cases, however, we did not find a sensitive dependence on the choice of 
potential. 
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Figure 5. Evolution of the variable x oc cb in the decaying dark energy model using the SUGRA 
potential, x passes through zero and becomes negative at In a = —12.336 causing an instability in the 
power spectra. This behaviour occurs for interaction strengths larger than ~ 10 _70 mpi. 

7.1 Instability in the dark energy decay model 

Our computation shows that the decaying dark energy model was found to suffer from an 
instability in the perturbations. From Equation 3.13, we see that significant values of in- 
teraction strength act to drive the value of x towards 0, causing dx/dN and to become 
infinite. Below the critical value of interaction strength associated with this instability, there 
is no significant effect on the power spectra or ISW effect. An example of this instability is 
shown in Figure 5 in which the interaction strength A is chosen to be just above the crit- 
ical value. Starting at z = 10 8 , we see that x decreases to at In a = —12.3 where the 
instability develops. By considering very early times, we can effectively set the limit on the 
interaction strength to and rule out this model of dark energy. This is because even small 
interaction strengths would, at early times, cause the small amount of scalar field kinetic 
energy to quickly dissipate, leading to this instability. In general, regardless of the choice of 
quintessence potential, we find that the susceptibility of the dark energy decay model to this 
instability makes it unviable as a model for interacting dark energy. 

7.2 Dark matter decay 

The primary effect of the dark matter decay model is to alter the time at which dark energy 
comes to dominate the Universe. This is because dark matter decay is most effective at early 
times when matter dominates over dark energy. This means that the background evolution 
of our 'best-fit' models is sensitive to the interaction strength. For instance, an interaction 
strength A = 10 _62 mpi changes the dark energy density at z = by roughly 1%, but leaves 
the matter power spectrum and ISW effect virtually unchanged. 

Figure 6 shows the evolution of the x variable for interaction strengths in the range 
10 -63 — 10 mpi. We see that the dark energy density is sharply boosted at early times, as 
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Figure 6. Evolution of the energy variable x oc </> in the decaying dark matter model using the RP 
potential. We see that increasing the interaction strength A boosts the value of x at early times which 
can lead to an increase in late time dark energy densities. 

the energy is immediately transferred from dark matter to dark energy. In Figure 7, we show 
the evolution of over the same range of interaction strengths. We see the corresponding 
sharp rise in towards 1 at early times. This has the effect of redshifting away the excess 
dark energy at a rate pj, oc a~ 6 . This rapid redshifting means that for a wide range of 
V (4>), dark energy does not come to dominate the universe too soon and so this model of 
interactions can indeed be reconciled with observations. 

A number of previous investigations of this interaction have used simple parameteri- 
sations of w^. However, we see from Figure 6 that the complex redshifting behaviour of 
dark energy would be extremely difficult to reproduce with a simple analytic form of w^. 
In fact, an instability in this form of interaction was found in [14] when considering dark 
energy parametrized by = constant. We can see that in a w =constant model where 
w < 0, the early boost in the dark energy density will not be redshifted away relative to the 
dark matter energy density and dark energy will come to dominate the Universe too soon for 
any appreciable interaction strength. It is interesting that due to the redshifting behaviour 
described above, our quintessence models do not encounter such an instability. 

7.3 Scalar-tensor model 

The scalar-tensor model exhibits the most interesting behaviour of the three interaction 
models considered in this work. The sign of the interaction can be changed to represent the 
flow of energy from one dark component to the other since the interaction depends both on 
dark matter and dark energy. If we choose = Ap m <j) where A is negative, we encounter 
the same instability in the x variable which we found in the dark energy decay model. 
From Equation 3.13, we can easily deduce whether a model with oc x n would develop 
an instability or not. When n < 1, the interaction term becomes increasingly negative as x 
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Figure 7. Evolution of for the dark matter decay models presented in Figure 6. At early times 
Wfj) responds to the sharp increase in x by rising suddenly to a value of 1. See the text for further 
discussion. 

decreases, hence leading to an instability. For models where n > 1, however, the negative 
constant tends to as x tends to and so an instability is avoided. In this model n = 1 
and to avoid an instability we require that the sign of the interaction term represent flow of 
energy from dark matter to dark energy, i.e. = Ap m cj) with A > 0. This supports the 
conclusion in [27] which found no background instabilities in this model. 

We find that this form of interaction does not change the energy density of dark energy 
significantly at early times since oc cj> and <f> is small at early times. The energy transfer 
becomes significant only at late times when dark energy becomes dominant. As we increase 
the interaction strength we see an enhancement in the matter power spectrum on large scales, 
as shown in Figure 8. The power spectrum on smaller scales k > O.OlMpc -1 decreases by 
~ 6% relative to ACDM for the largest interaction strength shown. 

We have shown that a scalar-tensor interaction can lead to an enhancement in the 
matter power spectrum on large scales. In order to detect this intrinsic enhancement we 
must be able to distinguish it from other artificial and observational effects such as gauge 
artifacts, the scale dependence of the galaxy bias and redshift space distortions which can 
also cause enhancements to the matter power spectrum on large scales and mimic our result 
(see [51], [52], [53] and [54]). We will now address these effects in more detail and explain how 
to break the degeneracy between an intrinsic large-scale enhancement due to a scalar-tensor 
interaction and a similar enhancement caused by the effects above. 

The work by [51] showed that calculating matter power spectra using gauge-dependent 
quantities can introduce gauge modes which can cause artificial large-scale enhancement in 
the power spectrum. We have used gauge-invariant quantities in our calculations and so 
the large-scale enhancement we have found is not due to such a gauge artifact. The linear 
bias is a linear relation between fractional overdensities of galaxies and matter (5 g oc 5 m ) 
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Figure 8. The linear matter power spectra in the SUGRA quintessence model with dark sector 
interaction of the 'scalar-tensor' form = Ap m (f>. Large values of A cause a significant enhance- 
ment to the power spectrum on very large scales. The interaction also leads to a smaller ~ 6% 
difference between the power spectra of the largest interaction strength and no interaction on scales 
k > O.OlMpc" 1 . 



which is used to infer the total matter power spectrum from a galaxy survey. Work by [52] 
showed that the linear bias relation is scale-independent only in the comoving-synchronous 
gauge and using perturbations denned in a different gauge leads to the linear bias having a 
scale-dependence. They showed that if the linear bias is used as a scale-independent relation 
in a different gauge it can lead to an artificial large-scale enhancement in the matter power 
spectrum [52]. Failure to use the correct form of the linear bias relation in a particular gauge 
could artificially replicate the signal from a scalar-tensor interaction. 

Redshift space distortions are caused by the peculiar velocities of observed luminous 
matter. These small peculiar velocities produce an additional doppler shift in the observed 
photon frequency relative to the average redshift caused by the background expansion velocity 
of the region. For large-scale overdensities the surrounding infalling matter will have an 
additional component of their peculiar velocities directed towards the overdensity due to 
gravitational attraction. This leads to the region appearing compacted in redshift space, 
which causes an apparent enhancement to these large-scale overdensities (since redshift is 
used as a proxy for distance). 

Recent work by [51], [53] and [54] calculated the effects of redshift space distortions, 
magnification by gravitational lensing and distortions to the luminosity distances of sources 
on measured quantities such as the matter power spectrum. They calculate photon geodesies 
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Figure 9. The fractional difference in the ISW effect in an interacting model relative to no interaction 
using the SUGRA potential. Only the i = 1 dipole moment of the largest two interaction strengths are 
observable above the cosmic variance. The largest interaction strength leads to an order of magnitude 
enhancement to the CMB dipole moment relative to ACDM. 



in a perturbed universe using linear perturbation theory. These geodesies take into account 
the deflection and redshifting of observed photons through interactions with overdensities. 
They then use the perturbed geodesies of the observed photons to relate gauge-invariant ob- 
servables from galaxy surveys (like the matter power spectrum) to the cosmological quantities 
calculated in theoretical models such as ours. Using the results of [51], [52], [53] and [54] we 
can break the degeneracy between an intrinsic large-scale enhancement due to interactions 
in the dark sector and enhancements due to the artificial and observational effects discussed 
above. 

Recently [16] investigated the effect of an interaction of the form (6.10) on the back- 
ground evolution, matter power spectrum and halo mass function using quintessence models. 
They found that significant negative values of the interaction strength A could lead to an 
enhancement in the matter power spectrum relative to ACDM, in agreement with our find- 
ings. The authors did not use fully gauge-invariant perturbation variables or include the 
effects of momentum transfer in their investigation, although they did include the effect of 
baryons and non-linearities. They did not find the instability in the x variable when A < 0, 
corresponding to dark energy decaying to dark matter. We expect that this is due to the 
substantially different treatment of the perturbations and the quintessence potentials they 
used. 

The production of large-scale enhancement of the matter power spectrum would be a 
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very significant result if it could be resolved in future observations. Recent work by [55] used 
peculiar velocities to constrain the matter power spectrum and found that an enhancement 
to the matter power spectrum on large scales was favoured by the data. However, due to 
large errors, primarily from cosmic variance, the findings were also consistent with ACDM. 
Nevertheless, peculiar velocities could potentially be used to constrain the large scale matter 
power spectrum in the future when a greater number of accurate measurements are available. 

Finally, Figure 9 shows how the ISW signal changes with increasing interaction strength. 
We find that the interaction yields an enhancement in the ISW signal around only the lowest 
multipoles (since the dark energy spectrum peaks at small k corresponding to small £). Only 
the 1=1 moment of the two strongest interaction strengths plotted (A = 10~ 175 mpj and 
A = 10~ 1 " 8 mp 1 ) is observable above cosmic variance. The largest interaction strength leads 
to an order of magnitude enhancement to the CMB dipole moment relative to ACDM. 

The recent high redshift survey by [56] found evidence of an excess clustering on large 
scales (k < 0.01/iMpc -1 ) at 4cr significance from ACDM. Using photometric redshifts from a 
sample of SDSS galaxies they measured the angular power spectrum and found an excess in 
the lowest multipoles of the angular power spectrum. This is a promising result, since it can 
be interpreted as the effect of the dark sector interaction shown in Figures 8 and 9. After 
rigorous analysis of their data and biases they concluded that the large-scale enhancement 
was a real effect which hinted at an exotic form of dark energy. 

We find the recent evidence for enhanced power in P(k) on large scales to be encouraging 
in light of our findings. However, we can see from Figures 8 and 9 that even if there were 
significant increments in the matter power spectrum at these near-horizon scales, they are 
only detectable in the dipole moment of the CMB, a measurement which is dominated by 
the Doppler shift caused by our own peculiar velocity. Unfortunately, the enhancements to 
the matter power spectrum on scales k < O.OOlMpc -1 are well beyond the scales probed in 
current galaxy surveys (k ~ O.OlMpc -1 ), so we do not expect that a direct detection of this 
large scale enhancement via galaxy counts will be possible in the near future. 

Using the cross-correlation of the CMB and large-scale structure from galaxy surveys is 
likely to be one of the most effective observational probes of our results. This technique allows 
us to isolate the ISW effect from other contributions to the CMB at low multipole moments 
since perturbations at the surface of last scattering were small compared to the observed 
large scale structure today. Large scale structure is correlated with the ISW effect since the 
ISW effect is caused by the evolution of gravitational potentials. Using the cross-correlation 
allows us to measure the ISW effect and compare this to the enhancement predicted in the 
scalar-tensor model for large interaction strengths (Figure 9). The work by [57] calculated 
the cross-correlation using WMAP 7 and several galaxy surveys and found an excess ISW 
cross-correlation of la from ACDM. This excess could be explained by a scalar-tensor dark 
sector interaction. 

8 Conclusion and discussion 

In this paper, we have calculated the background evolution, matter and dark energy power 
spectra and ISW effect for three quintessence models with and without dark sector interac- 
tions. Our calculations are based on gauge- invariant perturbation theory originally developed 
to calculate inflationary perturbations. We feel that our calculations are more reliable than 
those obtained when dynamical dark energy is parameterised by phenomenological ansatz for 
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w, since in our approach the clustering of dark energy can be directly linked to the underlying 
scalar- field perturbations. 

We found that without interactions our quintessence models differ from ACDM by at 
most a few percent in observable quantities, meaning that any differences are essentially 
unobservable. 

We have shown how a covariant treatment of dark sector interactions can give rise to 
momentum exchanges in the dark sector, which many authors have overlooked. Our approach 
complements and extends the work of [14], who examined momentum exchanges in constant 
w models using the dark matter decay interaction. We demonstrate our techniques on three 
model of interactions, namely, the dark energy decay, dark matter decay and scalar-tensor 
type of interaction. 

The dark energy decay model was found to have an instability in the background and 
perturbation equations, in agreement with previous works on this interaction, and thus can 
be ruled out. In the decaying dark matter model, we found that its primary effect was to 
change the evolution of the background energy densities, with no observable differences in 
the power spectra or ISW effect. We explained why w =constant models, when associated 
with this interaction, will lead to an instability, which is absent when field dynamics are 
taken into account. 

The scalar-tensor model of interactions appears to be the most interesting and robust 
of all three. We showed that it is stable if energy flows from dark matter to dark energy. 
For this case, we found that the interaction enhances the clustering of dark energy and dark 
matter on very large scales (k < O.OOlMpc -1 ) and produced an enhanced ISW effect which 
exceeds cosmic variance only in the lowest multipoles. We also found that the interaction 
produces a significant enhancement in the matter power spectra on very large scales. 

We find the recent indications of enhanced large scale clustering via peculiar velocity 
measurements [55] and enhancement to the lowest multipoles of the angular power spectrum 
[56] encouraging, although, measurements on these scales are very difficult to make and are 
often dominated by cosmic variance and our own peculiar velocity. 

Overall, the prospects of constraining dark sector interactions via the matter power 
spectrum and ISW effect appear daunting. The enhancement to the matter power spectrum 
on scales k < O.OOlMpc -1 is too large to be probed by current galaxy surveys k ~ O.OlMpc -1 . 
There are certainly other observational probes, however, which may reveal the presence of 
interactions in the dark sector, including i) cross-correlation of the ISW with galaxy and 
quasar distributions [57-59], ii) the growth rates of large-scale structures at high redshift, 
[60-62], iii) effects of large super-horizon perturbations on the CMB [63, 64], iv) gravitational 
lensing of the CMB using EPIC [65]. We envisage that our calculation techniques can be 
adapted to explore these issues at least in the linear regime. 

Finally, there will almost certainly be interesting non-linear effects on the matter power 
spectrum and CMB anisotropies arising from dark sector interactions. However, analytic 
progress in this regime is extremely challenging (though not impossible [44]). iV-body sim- 
ulations such as those initiated by [43] may hold the key to understanding the non-linear 
effects of dark sector interactions. 
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